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How marine communities are affected by C0 2 -induced climate change depends on the ability of species to 
tolerate or adapt to the new conditions, and how the altered characteristics of species influence the outcomes 
of key processes, such as competition and predation. Our study examines how near future C0 2 levels may 
affect the interactions between two damselfish species known to compete for space, and the effects of 
declining habitat quality on these interactions. The two focal species differed in their tolerance to elevated 
C0 2 , with the species that is competitively dominant under present day conditions being most affected. 
Field experiments showed that elevated C0 2 (945 uatm) reversed the competitive outcome between the two 
species with mortal consequences, and this reversal was accentuated in degraded habitats. Understanding 
these complex interactions will be crucial to predicting the likely composition of future communities under 
ocean acidification and climate change. 

The impact of climate change on marine communities will depend on how species' differences in physiological 
tolerances and performance alter biological interactions within a changing resource landscape 13 . Rising 
atmospheric C0 2 over the last century has already led to a reduction in ocean pH (ocean acidification) and 
an increase in ocean temperature 4 . Dissolved C0 2 levels in the ocean are rising in line with atmospheric C0 2 5 and 
are projected to exceed 900 uatm from the current day 390 uatm by 2100 6 , while temperatures of the coastal 
ocean may increase by up to 3°C. Coral reefs are particularly susceptible to these changes as they have evolved 
under relatively stable environmental conditions and many organisms may live near their maximum physio- 
logical tolerance limits for temperature and C0 2 7,s . How reef communities are affected will be dependent on the 
ability of species to tolerate or adapt to the new conditions and how the altered characteristics of species influence 
the outcomes of key processes, such as competition and predation 9 . 

Elevated levels of C0 2 affect the acid-base balance of marine organisms 10 and alters the processing of sensory 
information for organisms as diverse as snails, crabs and fishes 1113 . For fishes this effect occurs through modi- 
fication to the functioning of neuroreceptors 14 , which changes the way individuals interact with their habitats and 
with predators, leading to increased levels of mortality, particularly at critical life history bottlenecks such as 
settlement 1516 . While research demonstrates major changes in behaviour and performance for isolated indivi- 
duals, it is unknown how these conditions are likely to influence the biological interactions that form the 
foundation of community dynamics. 

Ocean acidification and elevated water temperature also affect the live coral that is an important source of 
shelter and food for many coral reef fishes. High C0 2 leads to a reduction in seawater pH that decreases the 
availability of carbonate ions required for coral skeleton growth 17 . In addition, ocean warming causes corals to 
lose the symbiotic zooxanthellae that produce much of their required food 18 . This coral bleaching has been 
responsible for wide scale loss of coral, and combined with the effect of storms has resulted in widespread 
degradation of live coral and the fish communities that associate with them 19 21 . While shifts in fish community 
composition have been associated with coral degradation, it is unclear how degraded reef habitat will affect 
interactions among fishes that live in these environments in a high C0 2 world. 

Marine fishes display inter- and intra-specific variability in the extent to which they are affected by C0 2 -driven 
modifications to their environment. While many species require live bushy coral as a nursery habitat 19,22 , many 
have the flexibility to satisfy their resource requirements from a broader resource base and are not confined to one 
habitat (often classed as habitat generalists). Moreover, some fish species are more susceptible than others to the 
effects of elevated C0 2 16-23 , and even within species there can be high variability in the degree to which sensory 
systems are compromised 13 . This variability in the way species are affected modifies the selection profiles by 
predators on juvenile fishes under moderately elevated C0 2 (—800-900 uatm C0 2 ) 24 26 . One of the fundamental 
processes that affect the susceptibility of individuals to predation is how the prey within a multispecies assemblage 
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interact to modify individual risk. Interference competition, is par- 
ticularly important in affecting the small scale distribution of growth, 
body size, body condition and distribution in relation to predation 
risk 27,28 . While there are some studies demonstrating differences in 
sensitivity to C0 2 among damselfish species, it is unclear how these 
differences may affect behavioural interactions under natural condi- 
tions of predation risk. 

These effects will not occur in isolation, but rather as the coral 
habitats degrade fishes that inhabit these habitats will themselves be 
battling to adapt physiologically to the new resources and balance of 
senses available to them for risk assessment 2 ". Our study examined 
how C0 2 levels projected to occur in the next 70- 1 00 y may affect the 
interactions between two damselfish species known to compete for 
space at settlement 28 and how these interactions were modified in 
coral habitat degraded through bleaching. First, we reared wild- 
caught damselfish larvae of the two species (Pomacentrus amboinen- 
sis and P. moluccensis) in elevated C0 2 (945 uatm) or control water 
for 4 days, which we have previously shown to be sufficiently long to 
yield the lasting behavioural effects associated with C0 2 -induced 
neuroreceptor modification 1415 . Secondly, we transplanted treated 
fish to an array of natural reefs made of unbleached or bleached hard 
coral in a 2 X 2 design to test if exposure to elevated C0 2 and coral 
state interacted to alter the behaviour of fishes and their risk of 
mortality during recruitment to adult habitat. Because C0 2 impairs 
the senses used to detect and respond to habitat characteristics and 
assess predation 1530 , we predicted that the effects of C0 2 might inter- 
act with habitat state to affect behaviour and ultimately survival. Our 
results show that elevated C0 2 reverses the competitive outcome 
between the two species with mortal consequences and that this 
reversal is accentuated in degraded habitats. 

Results 

There was a significant difference in mortality of P. moluccensis 
among the four treatments (X 2 = 20.05, df = 3, p = 0.0002). Fish 
from elevated C0 2 on unbleached coral had approximately four 
times the survival of fish from the other three treatments (Fig. la), 
indicating that there was an interaction between C0 2 treatment and 
coral state. Mortality curves for P. amboinensis also differed from one 
another (X 2 = 36.7728, df = 3, p < 0.0001), with fish that had been 
exposed to elevated C0 2 having lower survival than those from cur- 
rent day C0 2 conditions (Fig. lb). Those exposed to elevated C0 2 on 
bleached coral died within 48 h, while over 50% of fish from present 
day C0 2 conditions survived 72 h (Fig. lb). 

The maximum distance fishes ventured away from the coral shel- 
ter was affected by coral state and differed between species, but this 
species-specific effect was largely reversed with elevated C0 2 (Species 
X State X C0 2 treatment: F 1)208 = 5.10, p = 0.025; Fig. 2a). On 
unbleached coral at ambient C0 2 , both species stayed a similar dis- 
tance from the coral. However, on bleached coral, P. moluccensis was 
over twice as far from shelter than P. amboinensis. The trends 
between species were largely reversed by elevated C0 2 conditions 
(Fig. 2b). On healthy coral P. moluccensis stayed closer to the coral, at 
a distance that did not differ from ambient C0 2 condition. In con- 
trast, P. amboinensis was located almost 4 times as far from 
unbleached coral than when exposed to ambient C0 2 conditions, 
and twice as far from the unbleached coral than P. moluccensis. 
These trends were accentuated when the patch reef comprised of 
bleached coral (Fig. 2b). The height of fish on patch reefs differed 
between species, but the magnitude of the effect changed with C0 2 
conditions experienced (Species X C0 2 treatment: F lj208 = 12.31, p 

< 0.0006). P. amboinensis occupied the lower part of the patch 
regardless of C0 2 treatment or coral state, while P. moluccensis occu- 
pied the highest part of the reef, but was significantly closer to the 
bottom when fish had been exposed to elevated C0 2 (Tukey's test, p 

< 0.05; Fig. 3). 



The levels of aggression exhibited by fish to their size-matched pair 
was strongly affected by coral state and differed between species, but 
this species-specific effect was reversed with elevated C0 2 (Species X 
State X C0 2 treatment: F 1208 = 11.757, p = 0.0007; Fig. 4). Under 
ambient C0 2 conditions, P. amboinensis was much more aggressive 
than P. moluccensis with the dominant P. amboinensis being more 
aggressive on bleached coral (Fig. 4a). In contrast, under elevated 
C0 2 conditions the pattern was reversed by species, with P. moluc- 
censis being most aggressive on bleached coral compared to 
unbleached coral, but always more aggressive than P. amboinensis 
(Fig. 4b). 

Discussion 

Elevated C0 2 affected the outcome of competitive interactions 
between species, with a complete reversal in the species that was 
the most aggressive and dominant competitor. This reversal in the 
outcome of interference competition was most dramatic on 
degraded, bleached coral. The mechanism underlying this change 
was the differential effect of C0 2 on the behaviour of the two species. 
A previous laboratory study has shown that P. moluccensis is less 
susceptible to exposure to elevated C0 2 conditions than its congener 
P. amboinensis 16 , who is the competitive dominant under present day 
conditions 28 . For fish exposed to either current day or future C0 2 
conditions, the competitive interactions indirectly had mortal con- 
sequences because the aggressive interactions pushed the subordin- 
ate into areas of higher risk away from shelter; a pattern accentuated 
on degraded coral. The behavioural changes found in our study 
suggest that elevated C0 2 may result in changes in the way species 
interact with each other and this will alter the processes that shape 
communities and ecosystems. 

Understanding how altered species interactions affect communit- 
ies and ecosystems is fundamental to predicting the impacts of cli- 
mate change and ocean acidification 1,31 . Previous studies on benthic 
communities have observed enhanced competitive dominance of 
algae over corals under elevated C0 2 3,9 . In contrast, our results show 
a reversal of competitive dominance, indicating that ocean acidifica- 
tion may lead to far more complex outcomes than would be predicted 
with just enhanced competitive effects. Competitive reversals 
between species should be commonplace under natural circum- 
stances 32 for organisms with strong competitive asymmetry to coex- 
ist, however they are seldom conclusively demonstrated. Steneck 
et al. (1991) found that depth determined which of two species of 
crustose coralline algae was competitively dominant, but dominance 
was reversed when grazing limpets reduce fitness of one species 33 . 
Under present day C0 2 conditions, McCormick and Weaver (2012) 
found that who won the contest between P. amboinensis and P. 
moluccensis was strongly asymmetrical and dependent on size, 
over-and-above other ecological factors such as prior residency 28 . 
The importance of size in dictating the outcome of competitive con- 
tests is a common finding for many organisms 34 . P. moluccensis 
juveniles only started to win contests when they were 1.2 mm larger 
than P. amboinensis. Because the winner of competitive interactions 
is based on body size, the stochasticity of recruitment at small spatial 
scales 35 and its pulsed nature within a tropical recruitment season 36 
will mean that neither species will always win or lose. This may lead 
to a destabilizing influence on the competitive dominance of any one 
species, making competitive exclusion less likely. 

It holds therefore that changes in the competitive outcome will 
also occur when the species change in their ability to defend a 
limited resource. Studies that have examined the effect of C0 2 on 
the behaviour of solitary damselfish juveniles demonstrate that they 
display a willingness to take more risk, they fail to learn predators, 
and as a result die faster than control fish 15,24,37 . However, the ability 
to cope with elevated C0 2 varies among closely related species. 
Ferrari et al. (2011a) exposed 5 different species of damselfish from 
the genus Pomacentrus to damage-released chemical alarm cues 
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Figure 1 | Influence of competition on survival. Effects of interspecific competition on the proportional survival of (a) Pomacentrus moluccensis and 
(b) P. amboinensis on unbleached or bleached coral habitats when exposed to elevated or current day C0 2 levels. Kapalin-Meier survival trajectories 
illustrate the different survival trajectories on patches of live unbleached or live-bleached Pocillopora damicornis hard coral. Individual fish were placed in 
size-matched pairs on the habitats (n = 25-29). Photographs: Mark McCormick. 



from conspecifics 16 . Fish held at current day C0 2 levels all 
responded with a fright response to the alarm cues, while those 
acclimated to higher C0 2 levels exhibited different species-specific 
responses from no response to the alarm cues, through to a fright 
response that was little different from the controls. Further, in a 
mesocosm study, Ferrari et al. (2011b) found that the relative mor- 
tality through predation of 5 species of damselfishes changed when 
both predator and prey had been exposed to elevated C0 2 24 . 
However, the lack of behavioural information meant that the mech- 
anism underlying the result was unclear. Our study shows that P. 
amboinensis exposed to elevated C0 2 also becomes less aggressive 
and avoids similar sized competitors that they would normally have 
dominated. This lowering of aggression, which underlies interfer- 
ence competition, is critically important for the fish to maintain its 
status in the plankton feeding community and will have repercus- 
sions for such ecologically important activities as obtaining food of 
the highest possible quality and the acquisition and defence of 
spawning sites. 

Species are likely to acclimate to some extent to the modified 
chemical conditions found in future C0 2 conditions. However, pre- 
vious studies have shown that the behavioural responses of juvenile 



damselfish are the same in individuals kept for 4 days at elevated C0 2 
and those reared from hatching at the same C0 2 levels 15,30 . Therefore, 
there seems to be limited prospects for within-generation accli- 
mation of behavioural changes caused by exposure to high C0 2 , at 
least during the critical early life history stages investigated here. In 
one of the only transgenerational studies to date, Miller et al. (2012) 
found that the routine metabolic rate, growth and survival of Fl 
offspring of the clownfish, Amphiprion melanopus, exposed to 
1032 (iatm returned to levels found in control fish as a result of 
non-genetic parental effects 38 . The generality of this result is 
unknown but likely to vary among species. Interestingly, a recent 
study has demonstrated that damselfish may have limited ability to 
compensate for C0 2 -induced behavioural changes through parental 
effects over one generation 39 . The ability of the species used in the 
current study to acclimate over multiple generations by non-genetic 
parental effects is currently unknown, however, it is likely that com- 
peting species who have been differentially affected by C0 2 will be a 
realistic situation in the future. 

C0 2 appears to affect early life history traits such as respiration 
and growth in a different way to behaviour. Currently it is unknown 
for any species whether behaviour is likely to adapt over generations 
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Figure 2 | Influence of exposure to elevated CO2 and habitat type on shelter use. The mean maximum distance that fishes ventured from the shelter of 
unbleached or bleached hard coral within 3 min behavioural observations for Pomacentrus moluccensis (white) and P. amboinensis (grey) that had been 
reared under ambient or elevated C0 2 conditions. Letters above the bars represent Tukey's HSD groups. Error bars are standard errors (n = 26-29). 



to elevated C0 2 . Recent studies of marine fishes and invertebrates 
have shown that at least part of the behavioural effect of elevated C0 2 
appears to occur through a change in the gradient of acid-base rel- 
evant ions across the neuron membrane 40 . This leads to a change in 
the operation of neuroreceptors (GABA-A receptors), switching 
them from an inhibitory to excitatory role. Research is required to 
determine whether this neural mechanism can adapt over multiple 
generations before the full ramifications of the present study can be 
determined. 

Behavioural effects of high C0 2 will not occur in isolation; ocean 
warming that affects the health of corals that fish use as shelter 20 also 
has direct effects on the performance of interacting fish species. A 
number of recent studies have found that juvenile fish that use live 
coral as a nursery habitat move further away from the shelter when 
the coral is bleached, and this exposes them to higher predation 
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Figure 3 | Influence of exposure to elevated CO2 and fish species on space 
use. Mean proportional height above the reef base within 3 min 
behavioural observations for Pomacentrus moluccensis (white) and P. 
amboinensis (grey) that had been reared under ambient or elevated C0 2 
conditions. Letters above the bars represent Tukey's HSD groups. Error 
bars are standard errors (n = 53-55). 



pressure in the laboratory 41 and field 28 ' 42,43 . The movement away from 
bleached corals was also found in the current study, with the sub- 
ordinate species being furthest away from shelter. This was in part 
driven by an increase in the aggression by the dominant species, the 
identity of which changed in response to elevated C0 2 . Ocean warm- 
ing may be expected to accentuate this effect by elevating routine 
metabolic rates 40 , increasing activity 44,45 , and the daily food require- 
ments 46 , though some of this effect may be reduced through devel- 
opmental and parental acclimation 47,48 . 

Our results represent a first step towards an understanding of how 
community processes may change in response to elevated C0 2 . 
Future studies will also need to incorporate C0 2 effects on predators, 
who themselves could be variously affected 49 . As demonstrated in the 
present study, competition can affect survival by altering an indivi- 
dual's exposure to predators, either directly through aggression 
pushing subordinates into riskier habitats, or indirectly by elevating 
hunger in subordinates and in so doing influencing their balance 
between foraging and vigilance 50 . Because both intra- and inter-spe- 
cific competition is critically important in influencing prey vulner- 
ability to predators investigation of the three-way interaction 
between dominant, subordinate and predator under elevated C0 2 
will be the next crucial step in understanding the regulation of com- 
munities in a high C0 2 ocean. 

Methods 

Newly settled ambon damselfish, Pomacentrus amboinensis (Pomacentridae), and 
lemon damselfish, P. moluccensis compete for shelter in the same live coral living 
space when they first settle from the plankton 51 . Both are planktivores and show 
similar preferences for live, healthy coral at settlement 28 . Interaction with the other 
species at the reef edge reduces the growth of both species compared to locations 
where the other species is absent 28 . Within 2 mo on shallow reefs the two species 
display a disjunct distribution patterns with P. amboinensis associated with the base of 
the reef and P. moluccensis with live coral at the top of the reef 28 . 

Settlement-stage larvae of P. amboinensis and P. moluccensis were collected over- 
night using light traps (for design see Meekan et al. 2001 52 ) moored in open water 
around Lizard Island (14'40° S, 145'28° E), in the northern Great Barrier Reef, 
Australia. Fish were sorted to species and transferred to 35 L aquaria supplied with a 
continuous flow of either control (present day C0 2 ; 410 Liatm) or elevated-C0 2 
seawater (945 ua.tm; see Supplementary file) for 4 days. This period of time has been 
found to be sufficient to elicit the full extent of behavioural effects of high C0 2 15 . Fish 
were fed 4 times daily ad libitum with newly hatched Artemia sp. 

On the day of the field trial, fishes from the control and elevated -C0 2 treatments 
were placed into individually labelled plastic bags of seawater and measured with 
digital callipers. For each C0 2 treatment, interspecific pairs of P. amboinensis and 
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Figure 4 | Influence of exposure to elevated C0 2 and habitat type on aggression. Level of aggression displayed by juvenile Pomacentrus moluccensis 
(white) and P. amboinensis (grey) 45-60 min after release in size matched pairs onto a patch of habitat (healthy unbleached or bleached live Pocillopora 
damicornis). Fishes were either reared under ambient or elevated C0 2 conditions. Aggression Index is calculated as Displays + (3 X Bites) - Avoidances. 
Letters above the bars represent Tukey's HSD groups. Errors are standard errors, n — 26-29. 



P. moluccensis were created such that the size difference between fishes was no more 
than 0.3 mm SL. Fish length was standardised as it has been found to be important in 
determining the outcome of interaction between newly settled fishes and their sur- 
vival 53 ' 54 . There was no significant difference in size of the fish between species used in 
the experiment (mean ± se: P. amboinensis 11.4 ± 0.03; P. moluccensis 11.5 ± 0.03; 
Fi,2i4 — 3.3, p = 0.07). To reduce transport and handling stress, fish in bags were 
transported to the field site in a 60 L holding tank of seawater {to reduce temperature 
fluctuations) under subdued light conditions. 

Patch reefs used in the field experiment were composed of one of two states of the 
bushy hard coral, Pocillopora damicornis: live unbleached coral or bleached coral. 
Bleached coral was produced by collecting healthy coral and placing it in freshwater 
for 4 minutes 55 . Over the next 48 h stressed coral lost their zooxanthellae and became 
pale in colouration but still had their tissue intact. An array of patch reefs (—18 X 15 
X 18 cm) was constructed on a sandflat, arranged 4-5 m apart and 3-6 m away from 
continuous reef. 

Both damselfish species naturally settle on patch reefs near continuous reef. In this 
habitat juveniles are exposed to a diverse range of predators that use a variety of 
feeding modes from ambush (lizardfish Synodus dermatogenys and small cods 
Cephalopholis microprion) to pursuit (dottybacks Pseudochromis fuscus and wrasse 
Thalassoma lunare). These fishes can be observed preying on juveniles that venture 
too far from shelter. 

Experimental protocol. Length-matched pairs of P. amboinensis and P. moluccensis 
were placed onto unbleached or bleached coral patches. Patches were cleared of any 
other fishes or invertebrates using a hand net prior to release. Experimental fishes 
were released on natural reefs away from the study area. A small wire cage (—30 X 30 
X 30 cm, 12 mm mesh size) was placed over the patch to allow fish to acclimatise to 
their new surroundings while being protected from predators. Cages were removed 
40 to 60 min after release of the fish between 10:00 and 12:00 h. Fish presence was 
monitored 2 to 3 times per day (i.e. after the initial acclimatisation period, the evening 
after release and the following morning) for —72 h. The field experiment and 
monitoring of survival were conducted over a continuous lid period during a 
recruitment pulse in late October 2012. Water temperature during the census periods 
averaged 25.5°C. Our previous studies with tagged fish have shown that settlement 
stage damselfish do not migrate from these isolated reefs and that any loss is most 
parsimoniously attributed to predation (e.g. 4244 ). This research was conducted under 
animal ethics approval from JCU (Approval A1737) and a research permit from the 
Great Barrier Reef Marine Park Authority (G12-34811.1). 

Behavioural assessment. The behaviour of each fish placed on the patch reefs was 
monitored 45 to 60 min after placement on the reef following the protocol of 
McCormick (2009) 42 . Briefly, behaviour of the focal fish was assessed over a 3 min 
period by a scuba diver positioned 1.0 m away from the patch with the aid of a 
magnifying glass. Six aspects of activity and behaviour were assessed: a) maximum 
distance ventured from the habitat patch; b) height above substratum {categorized as 
% of the time spent within the bottom, middle or upper third of the patch); c) number 
of fin displays; d) the number of chases or bites; e) number of avoidance episodes in 
response to a conspecific approach. Two additional composite variables were created 
to summarize several of the variables recorded and enhance ecological interpretation. 
Relative height on the patch was summarized as a cumulative proportion of the time 
spent at varying heights over the 3 min observation period, with the top of the patch 
taken as height of 1, mid-patch a height of 0.5, and bottom a height of 0. For instance, 



when a fish spends 50% of its time at the top of the patch, 20% at the middle and 30% 
at the base the relative height value would be: (1 X 0.5) + (0.5 X 0.2) + (0 X 0.3) = 
0.6. An aggression index was also created by adding the number of displays to the 
product of three times the number of chases/bites and then subtracting the number of 
avoidance events 28 42 . Our other studies have shown that a 3 min observation period is 
sufficient to quantify behaviour because it is consistent in the short term { 1 0 min to 5 
days 44 ). The variable measured had been previously shown to be sensitive to 
competitive interactions between these two fish species 28,43 . 

Statistical analyses. Survival (up to 72 h) offish among the C0 2 treatment (ambient 
and elevated) by habitat {unbleached, bleached) combinations were compared using 
multiple-sample Survival Analysis, which uses a Cox's proportional hazard model 
(Statistica 12.0 56 ). This regression model is not based on any assumptions concerning 
the nature or shape of the underlying survival distribution. Survival curves of each 
C0 2 treatment and coral state were calculated and plotted using the Kaplan- Meier 
product-limit method. The Kaplan-Meier method is a non- parametric estimator of 
survival that incorporates incomplete (censored) observations, such as those cases 
where censuses had to be terminated on trials prior to their completion due to time 
limitations of a field trip. 

Three behavioural variables, maximum distance ventured, relative height and the 
aggression index, were compared between species, coral habitat and C0 2 treatment 
combinations with three-factor ANOVA (Type III SS) followed by Tukey's HSD tests. 
Residual analysis was used to examine assumptions of ANOVA. 
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